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1 Introduction 

The study of particle production is still an active research, especially in early universe 
inflationary cosmology. The non-linear dynamics of quantum fields in fixed cosmological 
backgrounds reveal new phenomena when compared to a Minkowski analysis. Such an 
analysis allows to obtain the metric dynamically from the quantum fields propagating 
in that metric. The dynamics of the scale factor a{t) is driven by semi classical Einstein 
equations 

'^M. + ^^ ^M. + (higher curvature) = - (t ) . (1) 



SttG^ ""^ SnG^ 



Here G^, A^ are the normalized values of Newton's constant and the cosmological con- 
stant respectively and G ^^ is the Einstein tensor. The higher curvature terms are included 
to absorb ultraviolet divergences. Recently [1] we solved equation (1) numerically taking 
in the right hand side of Eq.(l) the contributions due to one loop quantum corrections 



plus the energy density due to particle production. The energy density of the produced 
particles was calculated using the method of the complex time WKB approximation 
(CWKB) and the energy density so calculated through particle-antiparticle rotation 
turns out to be p^/a^-, as if it acts like a radiation density term. The result is quite 
surprising. However, the numerical solution of (1) reveals that the universe evolves in a 
de Sitter phase, avoids singularity and at late time behaves like a radiation dominated 
or matter dominated universe depending upon the mode of particle production. It is 
now believed that the parametric and spinodal instabilities give rise massive particle 
production and is a result of coherent oscillations of inflaton as it evolves in time. The 
equilibrium dynamics or out of equilibrium dynamics, whichever one we prefer to use in 
(1), we need solutions of Dirac equation (considering only spin-^ fields) and solve numer- 
ically equation (1). Such kind of numerical solutions considering scalar field has been 
carried out by Anderson [2-5], evaluating (T) for massless and massive fields with a 
chosen vacuum for quantum fields. For out of equilibrium situation, one evaluates (T^^) 
for any state |) such that the quantum field operator $(a;) is written as 

<l>(x) = 0(a;)+^(x), (2) 

where iIj{x) is a new quantum operator. Here 

0(a;) = ($(x)) , (3) 

mx)) = 0. (4) 

For a given Lagrangian using Hartree or large A^ approximation one evaluates 0(x) 
and iIj{x) and calculates e and P using renormalized equations of motion for dynamical 
evolution (for details please see [6-12]. We intend to apply this technique considering 
$(a;) as spinor field. Leaving aside the fluctuation term iIj{x) in (2), we study in this 
work the evolution of the fields on a fixed FRW background. Before going into the text of 
the paper we note that the numerical solutions of (1), with massless conformaly coupled 
scalar fields contributing to < T^^ > through one loop quantum corrections [see ref. 1-5] 
supposed to be valid also in spinor case when curvature is much greater than the mass of 
the field, reveal that the universe evolves in a deSitter phase, bounces and then emerges 
into classical universe, depending upon the matter content to be fixed by the mechanism 
of particle production in the model. In the literature, though most works concentrate on 
scalar fields, the effect of spinor particle production contributing to < T^j, > are not well 
tackled. In our view at a very early stage consideration of spinor particle production 



would be important if one has to investigate the preheating and reheating mechanism of 
early universe, where reheating (a process that starts when inflation ends) temperature is 
solely determined by the mechanism of particle production. The numerical investigation 
is no doubt a way to understand and study this aspect but it is a hard nut to crack. 
The analytical treatment allows one to channel the numerical work and this is what we 
intend to take for our discussion. This is also a reason to keep A^ term in equation (1). 
For further details the reader is referred to [2-5], and for numerical solutions backed by 
analytical input see ref. [1]. 

Our aim is now to study the evolution with scale factor a(t) = (t^)"". The cases 
n = |, n = I and n = 1 correspond respectively to radiation, matter dominated and 
Schrodinger model backgrounds. The calculations have been carried out in all these 
backgrounds. We report here the calculations of n = 1 case and cite also the results 
of de Sitter background for comparison. This work is also motivated to find an answer 
to the emergence of p^/a^ term as mentioned earlier. We do not intend here to the 
solution of dynamical back reaction equation (1), rather concentrate on the role of fields 
on particle production scenario according to CWKB [13-18]. 

Let us now understand the basis of the CWKB as follows. We write down the tempo- 
ral equation for Dirac particle in FRW spacetime in Klein-Gordon form. After separating 
out the spacetime part we determine the turning points as usual from the resulting one- 
dimensional Schrodinger-type equation in time variables. The turning points are found 
to be complex in the conformal time. The turning points are identified as points where 
particle turns back, not in space but in time, which according to Feynmann-Stuckleberg 
prescription must be identified as an event of pair production. Now there is a mathemat- 
ical analogy between the Dirac equation in a time varying electric field and the scattering 
off a potential which varies in space but not in time. The unitarity relation obtained 
from charge conservation relates the refiection and transmission coefficient R^ and T'^ 
which refer to the refiection and transmission coefficient not in space but in time. Pair 
production is considered as the refiection of a positron from the vacuum and R^ is the 
pair production amplitude. It is also found that in the CWKB, the particle production is 
basically a resonance particle production and the essence of the phenomenon of resonance 
is the rotation of a free-particle solution to a free-antiparticle solution as if the Dirac 
current shows a smooth interpolation from — J to + J. The problem in the interpolation 
is how to choose the initial conditions i.e., the particle and antiparticle states. To have 
a clear understanding on how to choose the initial conditions of quantum fields, we turn 



to the explicit solution of Dirac equation in fixed background and construct the current 
J^ = ip'-f'^ilj to look into the turnings of the current (if it exists) with the evolution of 
time. This is the main content of the present paper. We find marked differences for the 
cases n = 1, |, I when compared to the de Sitter evolution a(t) = exp{Ht). 

Though the method of CWKB was originally proposed by us in the context of the 
behaviour of quantum field in curved spacetime, recently many workers [19-21] are study- 
ing this aspect. Recently [22] we have established the rotation of currents that we are 
advocating during particle production and find that the CWKB results are gauge invari- 
ant result akin to Schwinger's gauge invariant result related to particle production. The 
present work is an attempt to study this 'rotation of current' in the context of curved 
spacetime, in order to understand the interplay of Minkowski and Rindler vacuum in 
relation to particle production. 

The plan of the paper is organized as follows. In section 2 we discuss the Dirac 
equation in expanding spacetime. In section 3 we describes the basics of CWKB. In 
section 4 we discuss the solutions for (i) a{t) = a^t and (ii) a(t) = exp{Ht) along with 
the corresponding currents. In section 5, we discuss the numerical results. The de 
Sitter case is discussed with a view to understand the emergence of radiation dominated 
universe from the act of particle production in section 6. We end up with a concluding 
section. 

2 The Dirac Equation in Expanding Non-Flat 
Spacetime 

In curved spacetime the Dirac equation is taken to be 

[h''{x)df, - i7^(x)r^(x)] ^ = m^. (5) 

where 7'^(x) are the curvature dependent Dirac matrices, r^(x) are the spin connections 
and m is the mass of the particle. 

The metric, we shall consider here, the Robertson- Walker type 



ds^ = dr - a\t) 



1 — kr"^ 



+ r^dO' + r' sin^ 9 dip' 



(6) 



that describes a homogeneous isotropic non-flat universe, where k = 0, +1, —1, corre- 
sponds to Euclidean space, a spherical space and a pseudo-spherical space respectively. 



From Eq.(6) we may write the metric tensor as 
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and its inverse will be of the form 
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(7) 



with p2 



kr"^. 



The afiine connections F" with a , n , u = 0, 1, 2, 3 calculated using Eqs.(7) 

(8) are as follows: 
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Eq.(6) may also be written as 



QjS — ±J ,,-Lj rvu (^^ GtX 



where the vierbein-components are 
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(12) 
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We may find out the relations between curvature-dependent Dirac matrices 7^(ic) and 
curvature-independent Dirac matrices 7^, having relations 7 7^ +7^7 = '^9 ^ and 7^7;,+ 
Itla = '^Vab respectively with /x, z/ = 0, 1,2,3 and r^^^ = 1, rj^^ = -1, r]^^ = for ? ^ j, 
with the help of the vierbeins L" by 



It^i^) = Kla, 



(15) 



and henceforth the curvature-dependent Dirac matrices 7^(0:) may be calculated from 
the relations 

l^{x)=g^^''j,{x). (16) 

Two sets of curvature-dependent Dirac matrices through the Eqs. (15) and (16) 
obtained in terms of curvature-independent Dirac matrices are as follows: 

7o(a;)=7o> 7i(a^) = -7i, 72(a;) = «^72, 73(2^) = «^ sin^73, (17) 



and 
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7°(a;) = 7o, l^i^) = 7i, l^{x) = 72> l^i^) 

a ar 



ar sin 9 



73- 



The spin connections are given by the equation 

1 



r, 



7''(a;),7^(x). 



(19) 



where, the covariant derivative is denoted by "semicolon (;)". Also using Eqs.(9), (10), 
(11), (12), (17) and (18) we find the spin connections through Eq.(19) as 



ro = o. 



1 d 
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1 2 p^«^i' 



{ar%% - P7i72) , 



Ts = 2 («^7o73 sin 9 - p7,73 sin d - 7373 cos Q) . 

Thus the Dirac equation (5) takes the form through the Eqs.(18) and (20) as 
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Using the unitary transformation S such that 



m 



^ = 0. 



where 



and its inverse 



^ = 5^-, 



-5=2 {%%+%% + %% + 1) e-5^^3^1 e-2^^i72, 



(20) 



(21) 

(22) 
(23) 



S-'=e-^^''^^^e'^'^'^'^^^{l-^,% -%%-%%), (24) 

the Dirac equation finally will be of the form 

dt + l-+ im^, - - JTf .^ + (p - l)ar (dr + i) 11 ^ = 0, (25) 

2a a I V r/ J J 

where, we have used 7. = 7oQ;j, ai = 7o7j and the transformation of the curvature- 
independent Dirac matrices in spherical polar coordinates by 




sin^ cos(/9 sin^ simp cos 6 
cos 6 cosip cos 6 sin(/9 — sin^ 
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(26) 



3 Basic Principles of CWKB 



Let us consider a one dimensional Schroedinger equation, not in space but in time. 



-tP + uj\t)tlj = 0. 



(27) 
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In CWKB we consider t to be a complex variable and assume uj{t) has complex turning 
points given by 

Lu\h,2) = 0. (28) 

Defining 

S{tf,U)= f' oo{t)dt, (29) 

the solution of (27) in CWKB is written as 

i){t) ' exp [iS{t, to)] + R exp [-iS{t, to)] • (30) 

Here to and t are real where to is arbitrary and to > t. In (30) R is given by 

-iexp[2ig(ti,to)] 

l + exp[2^S(ti,t2)]' ^ ' 

where ti,2 are the complex turning points determined from (28). The interpretation of 
(30) and (31) is as follows. In (30) the first term is the direct ray. It starts from to > ti, 
moving leftward arrives at a real t < to. The second term in (30) corresponds to reflected 
part. A wave starting from to reaches the complex turning point ti and after bouncing 
back from ti reaches t. It is represented as 

(-2)exp[25(ti,to)-25(t,ti)]. (32) 

The contribution (32) is then multiplied by the repeated reflections between ti and t2 
and the multiple reflection is written as 

OO -1 

^ \-z exp {tS{h, t,)}f^ = — (33) 

^0 l + exp[22^(ti,t2)] 

The combined contributions (32) and (33) comprise the second term of (30). For conve- 
nience we have neglected the WKB pre-exponential factor throughout. 

Let us now obtain a dynamical picture of particle production from the foregoing 
discussion of reflection in time. Consider a potential V ~ exp {—iuot) to supply energy 
for the pair (e+e^) creation. In Feynmann's space-time diagram, we represent it as in 
Fig. 1(a). 

\/^(e+,Ee+)T+(e-,^e-)T- (34) 

Here the positron of energy Ee+ and the electron of energy E^- both are moving forward 
in time out of the potential site and is represented by the arrow beside them. According 



to Feyninann-Stuckleberg (F-S) prescription, the negative energy particle solution prop- 
agating backward in time as being equivalent to positive energy anti-particle solution 
moves forward in time. Using this prescription, (see fig. 1(b)) 

F.S 

(e-,Ee-)T =^ (e+,-Ee-)i 

prescription 

i.e., a negative energy positron {Ee+ = —Eg-) moving backward in time (represented by 
down arrow) meets the potential site V (that acts as turning points) and moves forward 
in time. Thus F-S prescription applying to (34) now reads 

F.S 

(e+,Ee+) T +{e-,E,-) T =^ (e+,i?e+) T +(e+,-i5e-) i • 

prescription 

Thus pair production can be viewed as a process of refiection in time. For further details 
the reader is referred to [18] and also [13] and [14] to find the equivalence, 

[pair production amplitude | = [Reflection amplitude [. 

If we now consider exp [-\-iS{t, to)] as antiparticle solution, the reflected component is 
interpreted as a particle moving forward in time. This is the Klein paradox-like situation 
not in space but in time and R is interpreted as pair-production amplitude. The essence 
of (30) is that there is no particle at t ^ — oo,i.e., 

^m, _ exp[iS'(t,to)], t<0 (35) 

but at t ^ +00, (35) evolves into 

^»" + _L^ ^^P [^'^(^' ^o)] + ^ ^^P [~^'S'(t, to)] . 

We applied the results (30) and (31) in various expanding spacetime [12-14] with re- 
markable results. 



4 Solutions of Dirac Equation 

In case of flat spacetime i.e. k = O(orp = 1), we will, here, flnd out the solutions of the 
Dirac Eq.(25) and it will take its form as 



■ 3 a . 1 — > ^' 

Ot + -- + im% a . V 

2a a 



ij = 0. (36) 



Let 



^jk.x 
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/(k,t) 



{27rf'\g{k,t) 

and thus we get for two-component spinors /(k, t) and 5'(k, t) the coupled equations: 

3 a 



(37) 



and 



dt + -- + im] f k.ag = 0, 

2a / a 



dt + tm] g k.af = 0. 

2 a / a 



We get in turn the uncoupled equation for /(k, t) as 

/^n la 1 d^ . d 2 ^^\ 7 

d; h -^ + tm- + m^ + — Wi = 0, 



2 a 4a2 



where we have used 



f = -jh. 



(38) 
(39) 

(40) 
(41) 



4.1 Case-I: Schrodinger Model 



In Schrodinger model, we find 
Then Eq.(40) becomes 



a(t) = a^t. 



rfi , k a +- im 2 
* t t 



h = 0. 



(42) 



(43) 



After changing the variable from t to z = —Am t ot z = ±2imt, we get, from Eq.(43), 



dz^ 4 z z^ 



/i = 0, 



which is Whittaker differential equation having two independent solutions: 

h^{z)=W ,,(±22mt). 

Thus the upper components of the wave function in (37) become 



(44) 



(45) 



L 



:h^{z) 



2,2 , 

a^t ± 



W^,^,{±2imt) 



(46) 



10 



and the lower components will be 



g± = yi 



e [ K -K 



(47) 



where 
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and /Cj^ = /c^ ± ik^ . But Whittaker functions W^ satisfy the following identities: 
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We use the first identity for /t = +|, 2; = 2imt and for second identity k 
z = —2imt. Thus Eq.(48) gives, using Eqs.(49) and (50), 
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Denoting the constant spinor components of / by and 

solutions of the form (37) will be 
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the four independent 
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(55) 
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^ik.x 



i^. = N, 



(27r)3/2 a2t2 
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W, ^, {-2imt) 



(56) 



It is interesting to see the asymptotic form of these solutions for large times. Since 
^Kn{^) ~^ -2^6"^/^, for — ^ < arg2; < ^, we obtain the normalized solutions as follow: 
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Where we have determined the normalization constant A^. in such a way that asymp- 
totically, i.e., in the flat-space limit, we have the usual 5(k — k') normalization of the 
electron's wave function. The norm of ip is defined by 
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This procedure gives 
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(62) 
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4.2 Case-II: de Sitter spacetime 

In case of de Sitter spacetime, the function a{t) is given by 

a{t) = e^*. 
In a similar way, as in §4.1, we get four asymptotic normalized solutions as 
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We now evaluate the Dirac currents for these cases I and II. 



1 J.-i{z) 
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J J-u+l{z) 
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(67) 



Dirac Current in Case-I: 



Particle production in Schrodinger model may easily be interpreted from the nature of 
the curves of Dirac current in a given spacetime. Using Eq. (5) 



^l''{x)d^ - i-f''{x)T (x)] ip = m^. 



(6J 



in an external field A^ and its conjugate we can express the Gordon decomposition of 
current as follows [23]: 

1 



jt^ = ■^p■y^^■^p 



2m 






(69) 



13 



For k = i.e. in the flat spacetime the metric (6) reduces to 

ds"^ = dt^ — a?{t) {dx\ + dx\ + ^0:3) , 
and also the space-dependent Dirac matrices (18) get the following forms 

7°(a;) = 7o5 7^(a^) = — 7i, l'^{x) = — 7^, 7^(x) = — 73. 

a a a 

Using this equation (71) we have the current densities in components as 



^° = 2^"^^ (^^^■^«^) + i^ 



i dg +eAg 
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i a +e A, 



2m a^ 



(73) 



Now using equations (72) and (73), taking A = 0, we get the same results for current 
densities j^, jj, j^, and jg for the solutions tp-^ and ■^/'j given in Eqs.(57) and (58). These 
are as follows: 



J, = c,T + c'T , J, = -ex , j; = -C,T , J3 



-ar 



(74) 



and for the solutions ip.^ and ip^ of Eqs.(59) and (60), the components of the current 
densities get their forms as 



io ~ ^0-^ + ^0 ' ii ~ ^1^ ' ^2 ~ C^2-^ ? is ~ C^3-^ 5 



(75) 



where C^^ 
with T = 



c 



a 



a 



^vr^a? ' ° 327r3m2a5 ' ^ Sn^ma^ ' ^ 87r3ma 



^ and C, 

5 ■' 



^TT^ma^ 



t • 



Dirac Current in Case-II: 



Similarly the four components of the current densities, in de Sitter spacetime, will get 

their forms as 

Jo = C,T' 

j^ = C (k^ sinh ^-k^ sinT) T' 

j; = C (k^ sinh ^ + k^ sinT) T' 

3, = Ck,T%mh ^ 



(76) 
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or in brief, 

Jo = CJ' 

j,=CAl-DsmT)r 

j, = aAl + D,smT)T' 

Js = C,T' 
with the solution ^, and where, T = f e-^^*, Co = g^f^, C = ^^^^skr!tos(^i^) ^ ^^ = 
Ckismh{7rm/H), C2 = Ck2smh{7im/H), C3 = CA;3sinh(7rm/iJ), D = ^^ ^^,,^i^^^/H) and 

"■2 
Let us mention some arguments for citing the Gordon decomposition, though it is 

not required to obtain the expression of currents in the present work. We also calculated 

ip'-f^ip using the exact solutions of Dirac equation and find that this result obtained by 

us is coincident with the Gordon decomposition (73). In Barut's work [23] this Gordon 

decomposition is wrongly calculated and has forced us to check this result. The Gordon 

decomposition is carried out to investigate the effect of spin towards the rotation of 

currents which is extremely important in understanding the particle production leaving 

aside the complicacy of vacuum state definition in curved spacetime. 



5 Interpretation of Currents 

In order to interpret the currents qualitatively, we carry out the numerical analysis. For 
Schrodinger model as well as de Sitter spacetime we observe that j^ = i^'j^ip varies as 
■\. For the case-I, the nature of the curves for j^ , j^ , j^ are of the same form for the 
set of solutions ■?/'^ and ip^ . Similar results are also obtained for ip,^ and ip^. We observe 
that Jiiip-^ 2) — ~Jii'^3 4,)^ though j^ remains the same for both the sets. For convenience 
we take C^ = C^' = C^^ = C^ = C.^ = 1 and the curves are reproduced in figure- (2a) for 
Jo , figure- (2b) for the first set of solutions and figure- (2c) for the second set of solutions. 
The curves in figure- (2a) and figure- (2c) are of the same nature. 

The current in Schrodinger model has an interesting properties. From the figure- (2b) 
and figure- (2c) it is clear that the current of ip^ and ip^ is negative whereas for ip.^ and ip^ 
it is positive. So from equations (57) to (60) it is evident that the turning (also shown in 
figure- (2d)) (ip-^^jip^) -^ {ip^,ip^) would lead to particle production. Though the situation 
here is less attractive and interpretative compared to de Sitter spacetime, the plot of 
current density versus time provides interesting way to understand particle production 
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through particle-antiparticle rotation. Though we do not get a definite conclusions, we 
beheve that the particle production should somehow be dependent on initial states and 
the initial conditions would greatly effect particle production as well as the evolution 
of the universe. It is worthwhile to mention that the behaviour of current in radiation 
dominated universe is very much similar to Schrodinger model. It is evident from the 
curves that if there are some fluctuations that take the system out of equilibrium there 
will be oscillations between curves I and II in flg.-2(d). 

In case of de Sitter spacetime the forms of jo, Ji, J21 ^"^^ Jz ^"^^ of the same form as that 
of Eq.(77) when we consider any one of the solutions ip^, ip^ or ip ^ separately. Though 
there is some differences in constants, only we look into the behaviour of currents with 
time. The nature of the curves are of the same form when the current density j\, j,,, or j^ 
are plotted against time t . For convenience we take Cq = Ci = €2 = C^ = D = Di = 1, 
2k = H = 1 and the curve is reproduced in the flgure-(2e). The oscillating time- 
dependent current densities j\ and j^ always remain symmetrical about the exponen- 
tially decaying one j.^. The nature of the curve of j^ is of the same form as that of J3, 
not shown in the flgure. A single oscillation indicates particle production (though the 
nature of the vacuum is unspecifled). Repeated oscillations indicates particle produc- 
tion (at 2a, 2b, 2c, ...) and annihilation (at 2d, 2e, 2f, ...) in course of time. Without 
going into the details of quantum vacuum, it can be said that repeated oscillation is an 
indication of generation of photon along with particle production. The occurrence of 
oscillations at t large negative indicate the existence of repeated reflected trajectories 
as is demanded in our CWKB approach [13, 14,17,18]. The emergence of radiation dom- 
inated universe may thus be explained considering production and annihilation (trough 
and crest in the Fig.(2e)) as a mode for generation of photon at late time. In view of 
our work on de Sitter spacetime [15] the oscillating behaviour of de Sitter current is 
very interesting and worth-pointing in the sense that the emergence of radiation domi- 
nated universe might also also help understand the baryon-asymmetry problem in early 
universe. 

The numerical parameters adopted in this work are used only to have a quantitative 
understanding of photon production in the early universe. We will publish shortly the 
results taking parameters that corresponds to early universe situation. 
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6 Photon Dominance in de Sitter Spacetime 

Before we discuss the production of huge number of photons due to particle production 
in de Sitter spacetime, let us consider the arguments of Bernido [24] in this respect. The 
Green function for a relativistic scalar particle is expressed as 

G(x", X- m) = {thy' / exp( --) K{x", x'; A), (78) 

JO 2 

where x = (r, t) and A is a time-like variable. Antiparticles in the path integral approach 
are identified by particles which move backwards in ordinary time such that dt/dX < 0, 
with A always increasing. Pair production is understood as follows. Let A{r,t,X = 0), 
B{r', t,\ = \g) and C{r^,t^, A) be three spacetime points. Let a particle leave the point 
A and move backwards in time, scatter at C, turn back and move forward in time to 
arrive at B. Here r < r^ < r' and t^ < t. We note that here the path from A io B 
represents a particle moving backward in ordinary time t. In Feynmann's language we 
interpret it by saying that an observer at time t sees both a particle and antiparticle 
coming from C. 

In de Sitter spacetime, Eq.(78) has been evaluated in ref.[24]. It is found that the 
Green function yields no bound state spectrum and it exists provided t" > t' . It is to be 
noted that the spacetime points x'{= r',t') and x"{= r",t") can be viewed as the points 
where a particle and antiparticle can be detected. The consequence of the constraint 
t" > t is interpreted as follows. 

(i) For a given t, such that t' < t < t", an observer can only detect a particle and not 
its antiparticle. This gives rise to an observed matter universe. 

(ii) To obtain huge number of photons compared with the baryons as observed today, 
consider the particle starting at point x' = (r',t'). It may follow a path that jogs 
up and down (with respect to t coordinate) and a huge number of photons is thus 
generated from pair annihilation (i.e., creation at the trough and annihilation at 
the crest). 

The arguments placed in (i) and (ii) though seem encouraging but have some short- 
falls. What makes the particle jog up and down? A plausible answer comes from the 
CWKB approach where the mechanism described in (i) and (ii) fit very nicely. The 
details will be placed elsewhere and here we mention the essential aspects. The second 
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order Dirac equation or the equation (40) shows turning points at 

^^±,(,„^f), ,79) 

when we put a(t) = exp(ift). 

According to CWKB, the particle starts from t. and arrives at t with t < t. without 
any reflection. We call it direct rays. There is also a reflected part. A particle starts 
from t^ and moving backwards suffers a reflection at one of the turning points, say 
tj = +i (m + ^), and moving forward arrives at t. This is identified as pair production. 
The reflected part is enhanced by the repeated reflections between the turning points, 
in complex t-plane, given by equation (79). As like the arguments of Bernido [24], 
the particle jogs up and down resulting in a huge number of photons. The remarkable 
equivalence of CWKB with the standard techniques of particle production has already 
been established [15,16,18]. In CWKB, the repeated reflections gives rise to thermal 
spectrum and a reflection in time is considered as a rotation of currents i.e., —J to +J. 
In CWKB, we also find an answer to the Bernido's plausible assumptions. The matter- 
antimatter asymmetry can also be explained from our approach. However, we do not 
consider it here. Before ending, it is justified to ask whether the turning points given 
by (79) have anyway been reflected in the expression of Dirac currents in the curved 
spacetime. 

Consider the situation just after Big-Bang. We take k <^ Ha and m <^ H , the 
plausible conditions in inflationary cosmology. Under these conditions the expressions 
of the currents in Eq.(76) give 

- ~ ± -^ m, 
a \ki J 

quite consistent with Eq.(79), if we consider the case m'^ < (broken symmetry). The 

condition k <^ Ha correspond to long wavelength fluctuations greater than the Hubble 

radius, indicating the emergence of spinodal instabilities, a crude picture of which is 

reflected in the expression of Dirac current. However, the situation is not so simple 

and straightforward as mentioned. We need here to understand the interplay of vacuum 

energy and inflationary universe. The reader is referred to ref. [25] in this respect. 

It should be pointed out that in a realistic approach, huge photon production should 
follow from high temperature nonperturbative QED, but in curved spacetime a complete 
theory for which is still lacking. The work of Lotze [26] may be consulted for simultaneous 
creation of e'^e^ pairs and photons in RW universe. 



7 Conclusion 

The present work reveals that 

(i) the quantum vacuum at an early epoch starts from a state to be characterized 
at a given time, 

(ii) during the evolution of quantum vacuum there is rapid particle production 
and annihilation, of course in expanding spacetime. 

(iii) amplitude of oscillation decreases with time indicating a gradual attainment 
of stable ground state ( stable vacuum ), 

(iv) at late time particle production ends. 

In view of the above results it would thus be interesting to study the behaviour of 
Dirac currents in other expanding spacetimes. Rapid oscillation in the current is a char- 
acteristic feature of inflationary spacetime and may be introduced through an inflaton 
field or through a oscillating scale-factor (as is found in Starobinsky-type model without 
phase transition). We hope to discuss, in future, particle production, back reaction, 
avoidance of singularity [see 1], catastrophic particle-production through parametric res- 
onance within our framework. 
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Figure 2: A plot of the Current Density(j) versus Tinie(t) in arbitrary unit; (a) for j,,, in 
case of both sets of solutions, (b) for first set of solutions ipi, ■ip2, (c) for second set of 
solutions ips, ip^, (d) superposition of figures (b) and (c) in Schrodinger Model and (e) 
in case of de Sitter spacetime. 
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